Ideal quantum gas in expanding cavity: 
nature of non-adiabatic force 
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We consider a quantum gas of non-interacting particles confined in the expanding cavity, and 
investigate the nature of the non-adiabatic force which is generated from the gas and acts on the 
cavity wall. Firstly, with use of the time-dependent canonical transformation which transforms the 
expanding cavity to the non-expanding one, we can define the force operator. Secondly, applying the 
perturbative theory which works when the cavity wall begins to move at time origin, we find that the 
non-adiabatic force is quadratic in the wall velocity and thereby does not break the time-reversal 
symmetry, in contrast with the general belief. Finally, using an assembly of the transitionless 
quantum states, we obtain the nonadiabatic force exactly. The exact result justifies the vahdity 
of both the definition of force operator and the issue of the perturbative theory. The mysterious 
mechanism of nonadiabatic transition with use of transitionless quantum states is also explained. 
The study is done on both cases of the hard-wall and soft-wall confinement with the time-dependent 
confining length. 

PACS numbers: 05.20.Dd, 51.10.+y 



INTRODUCTION 



The nonadiabatic transition in the quantum gas of 
non-interacting particles confined in an expanding mi- 
croscopic cavity is the origin of the nonadiabatic force 
acting on the cavity wall. Let's consider non- interacting 
Fermi particles confined in a billiard with a moving wall. 
The wall receives the forces from Fermi particles in the 
billiard. Under the condition that whole system consist- 
ing of Fermi particles and a moving wall keeps the energy 
conservation, the work done on the wall by the force is 
supplied by the excess energy due to the energy loss of 
Fermi particles showing the non-adiabatic transition. In 
this way one can conceive both the adiabatic and nonadi- 
abatic forces. In the adiabatic limit, the adiabatic force 
due to the quantal gas on the cavity wall is proportional 
to the derivative of the confining energy with respect to 
the cavity size. What is a characteristic feature of the 
nonadiabatic force when the cavity wall is moving? The 
general belief is that this force should be linear in the wall 
velocity, breaking the time reversal symmetry. In fact, 
in compound systems like molecules where two kind of 
coordinates with different time scales coexist, the Born- 
Oppenheimer approximation leads to both the adiabatic 
and nonadiabatic forces acting on the degree of freedom 
characterized by the slow time scale, and the latter force 



is linear in the velocity of the slow degree of freedom and 
is called an irreversible or frictional force In the 

case of the expanding cavity, however, the Hilbert space 
as well as the domain of Hamiltonian is time-dependent 
because of the time-dependent length scale of the cavity 
confining particles, which requires a deeper insight into 
the nature of the nonadiabatic force. 

The investigation of the above subject was started by 
Hill and Wheeler in 1952 in the context of nuclear 
physics. The nature of the nonadiabatic force on the 
cavity wall were intensively studied by Blocki et al. [Hi- 
Wilkinson developed the extensive study on this subject 

01 in the context of energy diffusion and of random 
matrix theory assimilating the chaotic motion of parti- 
cles inside the cavity, which was followed by other groups 
[s- 11|. Most of these works regarded the force due to the 



quantal gas as conjugate to a time-dependent wall coor- 
dinate. However, the definition of the force operator is 
not clear at all in the case of the hard wall. In fact, one 
cannot define the force operator by using a given Hamil- 
tonian for the billiard with a moving boundary. Further, 
without the valid definition of the force operator, essen- 
tial results so far would be highly questionable. 

In this paper, we introduce the force operator with use 
of the time-dependent canonical transformation which 
transforms the expanding cavity to the non-expanding 
one. Applying the perturbative and exact theories, we 
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the wall receive the extra force depending on its velocity, 
which comes from the nonadiabatic transition occuring 
in the quantal gas. However, the definition of the force 
operator is far from obvious in the case of a moving hard 
wall. Below we shall define the force operator in two 
ways. 

A. Classical force and quantization 



FIG. 1: Moving wall confining the quantal gas. L{t) is the 
time-dependent size of the cavity, and F(t) stands for the sum 
of adiabatic and non-adiabatic forces. 

evaluate the non-adiabatic force whose nature thoroughly 
differs from the conventional frictional force. For com- 
parison we shall also investigate the case of the soft-wall 
confinement with the time dependent confining length. 

In Section |lTl we construct the force operator acting 
on the moving wall in an unambiguous way. We consider 
the case that the cavity wall suddenly begins to move at 
time origin (see Fig[T|) . In Section IIIIl within a frame- 
work of von Neumann equation for the density operator, 
we apply a perturbative theory to obtain the expecta- 
tion of the nonadiabatic force. In Secti on IIVI with use 

we ex- 
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of the transitionless basis functions 
actly evaluate the energy expectation to see the nature 
of the nonadiabatic force, and justify the validity of both 
the definition of force operator and the issue of the per- 
turbative theory. In Section |V] we investigate the case 
of soft-wall confinement by treating a tunable harmonic 
trap. Section IVll is devoted to summary and discussions. 
In Appendix |^ we investigate the exactly-solvable case 
when the expansion rate of the cavity is time-dependent, 
to see the universality of the assertion of the text. Ap- 
pendix [B] treats technical details of the integrals used in 
Section HVl 



II. 



FORCE OPERATOR 



1} "ri , gives rise to = 



which, with 

^ '22 

n n 



When a given cavity in 1 dimension has a size L and its 
wall is fixed, the force on the wall due to the quantal gas 
inside the cavity is defined by = 
use of eigenvalues En — 
And the contribution from all particles is expressed as 
F = J2Zi ^/(^n) where /(£;„) stands for the Fermi 
distribution function. At zero temperature, f{En) — 1 
for 1 < n < iV and f{En) = otherwise. This force 
is called the adiabatic force. When the wall will move. 



As a first step to provide the force operator, we show 
a kinetic evaluation of the force due to the classical ideal 
gas. Then the force is expressed as a dynamical quantity, 
and is thus quantized straightforwardly. 

Suppose that N mutually noninteracting particles with 
a common mass m are confined in a 3-dimensional (3- 
d) box whose 3 edges have a common length L. Along 
each of 3 coordinates, one wall is fixed at the origin, for 
instance, at x = 0, and another one is initially located, for 
instance, at x = L and begins to move with a constant 
velocity (i). The velocity of the wall L is assumed to 
be slow compared with the mean square velocity of the 
particles {vf). In the course of time evolution, the 
particles eventually become uniformly distributed in the 
box. The time t for a particle at a position x = Xi 
running to right with Vi{> 0) (to left with —Vi) to reach 
the moving wall is given hy t — and its average is 



t = 



Vi - L 



(1) 



The average time for a particle to come back to the initial 
position is given by 2t. On the other hand, the momen- 
tum change at each collision with the moving wall is given 
by 

Api = m{vi — L) — m{~Vi ~ L) = 2mvi. (2) 

Since the collision rate is given as the inverse of 2t, the 
force acting on the wall is given by 



N ■ 

yr-^ Vi - L 



^ 2L 

1=1 

( L 

I —mvi - —mvi 



(3) 



The first and second terms on the second line are the 
adiabatic and non-adiabatic forces, respectively. These 
forces are rewritten as 



^2 2 
^ rnvf_ ^ 

1=1 



(4) 
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and 



TV 



'^iTT^iLmv^+mViL) = -{^{xp + px)), (5) 



where p is the momentum. For each particle, the bracket 
evaluates the value of the position x and momentum p at 
the instance of collision. 

Then the non-adiabatic force is given as 



L , ^ 

Fnon-ad = - -jJ^KXp + px) . 



(6) 



As a dynamical quantity, Fnon-ad is invariant under the 
time reversal operation, since both the expansion rate L 
and momentum p change their signs. 

Let's quantize the non-adiabatic force obtained above. 
The force operator should satisfy the following condi- 
tions. 



a) In the classical limit, Fnon~ad agrees with 
-TT^ixp + px). 

b) Fnon-ad should bc Hcrmitian. 



c) Fnon-ad docs uot depend on the particle statistics 
(boson or fermion). 



Consequently the force operator should be 



non — ad ~" ^^J^ \ Xp -\- pX j . 



(7) 



The idea above is based on the phenomenological ar- 
gument with use of a 3-d box, but suggesting a promising 
expression of the force operator. Below we shall provide 
a rigorous definition of the force operator. 



B. Rigorous definition of force operator via 
time-dependent canonical transformation 

The original Hamiltonian H for the billiard with a 
time-dependent cavity size L{t) is given by (in unit of 
fi^/m ^ 1) 

1 ^2 



H 



2dx^' 



(8) 



We now see the expectation of H as given by 

imw, (9) 

where 1-0) is a solution of the time-dependent Schrodinger 
equation 



d 



(10) 



with a moving Dirichlet boundary condition: 

^{x = 0, t) = V(a; = L{t), t) = 0. (11) 

The expectation of the force acting on the wall is ob- 
tained by 

d 



F = - 



dL{t) 



{Mm)- 



(12) 



Noting ■§i\i^) = x^lV") = TKl^\'^^ ^^'^ Hermitian 
conjugate, Eq. reduces to 



- BM 

Hence the force operator is defined by 

OH 



F = - 



dL{t)' 



(13) 



(14) 



However, the original Hamiltonian H for the billiard with 
its time-dependent size L{t) does not formally include 
L{t) explicitly. Therefore there is no way to define the 
force operator directly by using Eg. (ITU) . 

To overcome this difficulty, we shall make the time- 
dependent canonical transformation of H related to the 
scale transformation of both the coordinate x and am- 
plitude of the wave function ?/;. This transformation, 
which was originally developed in the heat equation the- 



ory 
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131, is defined by [IJ-II6I 



i?i=e-^^(il-*;i^)e'^. 



(15) 



where 



U = -^{xp + px)\nL{t) = i (^x^ + hi L{t). 

(16) 

This canonical transformation leads to the scaling of the 
coordinate x, 



e-'^xe"^ 



xL{t), 



(17) 



where on the right-hand side the new variable x{= y) 
varies in the range < y < 1, which is time- independent! 
Similarly the amplitude of the wave function is scaled as 



4)ix,t) = e~'^^{x,t) = \/L'il;{xL,t), 



(18) 



so that the normalization factor of 4>{x, t) becomes L- 
independent. Finally the Schrodinger equation is trans- 
formed to 



ih— = Hi(, 
at 



(19) 
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with the new Hamiltonian 

d ihL 

4>{x,t) now satisfies the fixed Dirichlet boundary condi- 
tion 0(0, = = 0. Equation ^ with (gHl) is also 
available simply by replacing i/) and x and by-^0 and 
xL respectively in Eq- PH]) with Eq.®. 

Taking L derivative of Hi , we can rigorously define the 
force operator in the transformed space as 

1 ihL ( d 



F 



dHi 
1)L 



L^dx^ L2 V dx 



j^Pl - ^ i^Px + PxX) . 



(21) 



Now, carrying out the inverse canonical transformation 
{xL — 7> a;, etc.), we have the force operator expressed in 
the original space as 



L 



L dx"^ \ dx 



(xp + px), 



L 2L2' 
which certainly satisfies: 

(V-IFI^) = {4>\F[ 



(22) 



(23) 



The non-adiabatic term in Eq. (f22|) agrees exactly with 
the phenomenological result in Eq. (O. 

III. PERTURBATIVE THEORY OF 
NONADIABATIC FORCE 

In this Section we shall investigate the expectation of 
the force operator in Section|TT]perturbatively with use of 
von Neumann equation for the density operator and adi- 
abatic bases. The method is an extension of the Green- 
wood's linear response theory [2^. Let' s assume that 
the cavity wall is fixed with the cavity size Lq until the 
time origin t — and that it suddenly begins to move 
with constant velocity L at t > 0. 

In the equilibrium statistical mechanics, the expecta- 
tion of a given observable O is defined in energy-diagonal 
representation, as 

1 \ NT^ „ 1 



(O) = Tr O 



g/3(H-p) + l 



(24) 



In the near-equilibrium, the expectation value is eval- 
uated in terms of the density operator p 



P 



E 



a)uja{a\, 



(25) 



'^ith Sa'^a = 1, as 

(O) =Tr(p6) =^t^,0„ 



(26) 



We shall employ the original Hamiltonian H and coor- 
dinate X. The density operator p for the Fermi gas obeys 
von Neumann equation 



^n% = [H,pV 



(27) 



With use of adiabatic basis | n) , the matrix elements of p 
satisfies 

dp 

Pnm = ("1^1"^) + {h\p\m) -f {n\p\m), (28) 
at 

with 

{n\p\m) -f {n\p\rn) = ^ ((n|OP£m + P7ii{^\rn)) . (29) 

For the system under consideration, the instantaneous 
(adiabatic) eigenvalue problem is given by 

H{t)^l^n = EMt), (30) 

with adiabatic eigenstates and eigenvalues 



ibn (= {x\n)) — . / — -— sin ( — -— 
^" ^ ^ ' ''^ y L{t) \L{t) 

En — 



(31a) 
(31b) 



where we prescribed fi^/m = 1. 

Using Eq. pOI) . we can obtain the following formulas 



{£\m) = --7fm, 
{n\l) = {i\nr = -^jer., 



(32a) 
(32b) 



where 



7,™^ (-1)^+^+1^5^ (l-<5,„0. (33) 

Noting the pure-real nature of the adiabatic states in 
Eq.dHl]), we see {i\e) = {i\£) = 0. Hence we can put the 
diagonal element ^ee = in Eq. p3p . The von Neumann 
equation now becomes 



— (En — Em) Pnm ~ ( ^ UnPtm + ^ 



"ilmPni 



.l^n 



(34) 
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Eq. (l34l) can be solved perturbatively. Let' s assume 
the solution to be expanded in 0{L/L) as, 

P-m)^,M^9.(m]\-- (35) 



L(o) -'^{m 



with f{H) = -^^IjT^TJTTp^- 

One sees, for 0(1), 



fnm =0 {n = m). 



Therefore 



fn 



1 



g^j(£;„(o)-M) + 1 
Then, for 0{L/L), one sees 
En — E, 



(36) 



(37) 



gin 



ih 



+ 7nm/n), (38) 



where the result in Eq. (1571) was used. The solution of 
Eq. ([38|) is given by 



(l - e^^*) (/„ - U. (39) 



For a correction of 0{{L/ L)'^), the dominant contribution 
comes from the diagonal term satisfying 



g2nn — — ^ lin{giln + 51"^)- 
f. 

With use of Eq. 



(40) 



„ . „ „ N ^ . En - El 
gUn + 5ln<? = -^7en[Jn - Jl)— ^ Sm : 1. 



En -El h 



(41) 



So, using Eqs. pOI) and (HTI) . we obtain 



52« 



1-cos ^". ^' n. (42) 



Now let' s calculate the matrix elements of the force 
operator in Eq. ([2^ . Using the adiabatic bases, we find 

2 Z"^**) . frmTx\ ^ . /n7rx\ , 
Fmn — — r-r / sm — -— F sm — — ax 
L{t)J„ \L{t)J \L{t)J 



{n^f ^L{t) 



(43) 



Combining Eq. (gS]) with Eqs. ([371), dSH]) and (gH), the 
expectation value of the force operator becomes 

F={F)= Tr(pF) = Pr.mF^mn = Si + S2 + S3, 

(44) 



where 



Si 



(45) 



(46) 



^^(^(0))^ 
(L(0))2L3 

En — El _ 



t X 



(47) 



5*1 in Eq. (|45|) gives rise to the expression for the adi 
abatic force at finite temperature 

(n7r)2 1 



^ £3(0) e^(^"(0)-'') + 1' 



(48) 



while S2 in Eq. (gH]) and in Eq. (gT]) contribute to the 
nonadiabatic force. In Eq. (j48| . we have assumed the 
time range lies in 



h L 
AE«'«V 



(49) 



where the lower and upper limits of the inequality in 
Eq. ()49p imply the minimum resolution of time and the 
time necessary for the wall to move by order of L, respec- 
tively. This is a physically imposed assumption, which 
will also be employed below. 

S2 can be rewritten, using ^nmlmn = - (^iT^l^a 
and L'^{t){En - Em) = ^{n'^ - m^) and 
noting the fact that, for a symmetric func- 
tion s{n,m), Em^n^l*^'™) (1 ~ e-^'TTr^*^ 

2Em#«5('^'™)sin^ ( ^"27f'" ^) ^ Em#„s(":»T^). ^hcrc 
the final reduction is possible under the assumption in 
Eq. (gS]). As a result, 5*2 becomes 



^2- 



i6;i'£(o)^ ^ 



{fn fm) 



(50) 



The factor (/(£■„) - f{Em)) in Eq. ([50]) gives a con- 
straint under which the summation over m and n should 
be taken. Similarly S3 becomes 



^3 



32;i2(i.(o))^ 

7r2(L(0)) 



{In fm) 



(51) 
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To summarize, the nonadiabatic force is given by 



Fnon-ad — S2 + S3 — C 



with 



{fn fm) 



(52) 



(53) 



where C < since Fermi distribution function mono- 
tonically decreases with energy and fn — fm = f{En) — 
fiEm) < for n > m. Fnon-ad in Eq. is propor- 
tional to and does not break the time-reversal sym- 
metry, in marked contrast with the general conjecture so 
far. The result in Eqs. pS)) and (15^ will be confirmed 
by the exact analytical result in the next Section. 
Here we should give two comments: 

i) The first one is concerned with the level crossing. 
From an experimental viewpoint, we are considering a 
quasi-one-dimensional (1-d) hard- or soft-walled rectan- 
gular parallelepiped. In this case, the energy gaps be- 
tween sub-bands are large enough not to meet cross- 
ings among sub-bands. Therefore the dynamics within 
each sub-bands (e.g., the lowest sub-band) used in our 
scheme is guaranteed. As a more general case, one might 
consider a 3-d rectangular parallelepiped with the size 

one of whose walls is moving in x-direction. 
Then each adiabatic state is characterized by a set of 
quantum numbers (n^;, riy, n^) and the energy spectra as 
a function of might show level crossings among man- 
ifolds with different nx,ny and n^. If a confined particle 
is initially in a manifold with the fixed Uy and and 
the cavity expands only in x-direction, however, there 
occurs no transition among manifolds with different riy 
and Hz and thereby energy crossings do not affect the 
present dynamics at all. Finally, one can conceive an 
expanding 3-d spherical billiard, which has level cross- 
ings among manifolds with different angular momenta. 
Since there is no transition matrix element among differ- 
ent angular momentum states in the symmetry-keeping 
dynamics, however, the dynamics is free from the prob- 
lem of level crossings if a zero-angular momentum state 
will be chosen as an initial state, which again guarantees 
our scheme. 

ii) The second one is whether or not the expression 
for the force operator and the expectation of the non- 
adiabatic force quadratic in the rate of dilation under the 
"time-dependent" Dirichlet boundary condition (TDD) 
would be available directly from a variational method. 
Berry and Wilkinson (BW) [l3| considered the static 



(adiabatic) eigenvalue problem of a triangular billiard un- 
der the " time- independent" Dirichlet boundary condition 
(TID), to study the density of diabolical points and both 
shifts and splitting of level degeneracies. They had re- 
course to a degenerate perturbation theory with use of 
diabatic eigenstates at the degenerating point, indicating 
(in Appendix of their paper) that the off-diagonal energy 
matrix elements are zero for dilations. It is not easy to 
interpret the present dynamical result under TDD only 
in terms of the static one of BW under TID. The present 
work is concerned with a one-dimensional billiard with 
a time-dependent walls, where the adiabatic eigenvalues 
have no degeneracy and only the level shifts occur against 
the adiabatic dilation. In the context of the adiabatic 
(AF) and non-adiabatic forces (NAF) which are given 
on the second line in Eq.(|43|). the non-diagonal matrix 
elements of AF certainly vanish, consistent with BW, 
but matrix elements of NAF are new, whose counterpart 
cannot be found in the treatment of BW under TID. The 
force operator of BW is only concerned with AF. Non- 
vanishing matrix elements of NAF in Eq. (|43| are due to 
a dynamical contribution in Ea. (l22p coming from TDD. 
A mechanism of the absence of a term linear in the rate 
of dilation in the expectation value of the force operator 
in Eq. (j44p is not due to vanishing non-diagonal matrix 
elements of AF and can not be explained directly within 
the static framework under TID. It is caused by a sub- 
tle cancellation of the linear cross-coupling terms among 
the matrix elements of the force operator expressed as a 
series expansion w.r.t the rate of dilation and the density 
matrix expressed in the similar expansions in the frame- 
work of the extended Kubo-Greenwood formula. On the 



other hand. Berry and Klein 18| were once involved in 
the similar subject as the present one, but they showed 
neither the definition of NAF operator nor the expecta- 
tion value of the force as a power series in the rate of 
change of the scale size of the container. 



IV. EXACT ANALYSIS 

Exact solution of the Schrodinger equation with a mov- 
ing Dirichlet boundary condition due to the motion of a 
wall was found by Makowski et.al 19, 2^. The greatness 
of their work lies in that they discovered the transition- 
less basis functions where the adiabatic states are also 
the solution of the time-dependent Schrodinger equation, 
which recently received a great attention in the context 
of the shortcut to the adiabatic dynamics 2l|. With 
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use of their basis functions we can proceed to evaluate 
the nonadiabatic force exactly. After a brief summary of 
their results, we shall carry out this procedure. 

The system we are going to explore is described by the 
Schrodinger equation 

" "2 



ih- 



dt 2 dx^ 

where the wave function satisfies the moving Dirichlet 
boundary condition in Eq. pT|) . 

After the scaling of space coordinate x and wave func- 
tion il) by L{t) and ^yi^(t), respectively as 

X 



x^y = 



^ = VW)^^ (55) 
Schrodinger equation with the moving boundary becomes 

with a static Dirichlet boundary condition, i.e., 

^(0) = ^(1) = 0. 

The transformation above is nothing but the time- 
dependent canonical transformation described in Section 

im 

Then, applying the gauge transformation 

ih 

(t}{y,t) = G(j){y,t) = exp ' 



L{t)L{t)y'\ct>{y,t), (57) 



Eq. (|56p can be reduced to the Schrodinger equation for 
the time-dependent harmonic oscillator: 



(58) 



Eq. ([55)) can be solved exactly if the time-dependence 
of the boundary satisfies the following equation [l9| : 



L^L = const = --B^. 



(59) 



For a linearly expanding or contracting billiard with 
the constant wall velocity L = L{0), i.e.. 



L{t) = La + Lt, 



(60) 



the condition (15^1) is satisfied, B = 0. (A general case of 
B will be investigated in Appendix Ell 
The solution in this case is 



12 ihL ^ 



2h 



-r{t) 



I nnx 



X sm 



(61) 



where rit) is a new time variable defined by 

ds 



(62) 



/o L^s)- 

The solution (|5T|) is the transitionless state where the 
adiabatic state also serves as the solution of the time- 
(54) dependent Schrodinger ec 



mation, which recently received 
a renewed attention [21 1. An assembly of states in 
Eq. ([6T|) constitute the complete ortho- normal set. 

Let's obtain the adiabatic and nonadibatic forces act- 
ing on the moving wall which is confining the Fermi gas 
into the cavity, by evaluating the expectation of Hamilto- 
nian. Statistical weight factors (Fermi distribution) will 
be incorporated a posteriori. In the case of a linearly ex- 
panding cavity described by Eq. ([60)) . the initial state of 
a particle is given by ^^^^sin ) '^ith the eigenvalue 

El (0) = and the wall suddenly begins to move with 



constant velocity L. The solution of Eq. (IMI) can be ex- 
pressed in terms of the transitionless states in Eq. (pTj) 
as 



1 



ih _ 



^^Cn(t)fn{y), (fn{y) ^ V2sm{mTy), (63) 



where expansion coefficients are given by 

'-Tit) I (64) 



cl^\t) =cW(0)exp 



2h 



with 



ih _ 



cW(0) = 2 sin {Iny) sin {mTy)exv\^^ J L{0)L{0)y'jdy. 

(65) 

The average energy can be represented as 

rL{t) 

{E^'^\t))= il;*{x,t)Hi^{x,t)dx 
Jo 



1 

LJo 



riy,t)[-^]^P{y,t)dy 



1 hi , , h^P , , 

= + —lm{h) + -^h, (66) 

where Iq , Ii and I2 are respectively defined by 

lo^ [ I P dy 
Jo 

=27r^ ^ cl^)* {t)n' ^ c^^' {t)n / cos{n' Try) cos{mry)dy 
=n'j:\c(:H0)\'n\ (67) 
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-'0 



(68) 



hit) = f U P = Y^c^;}*it)Y^c\f\t).hin,n% 

(69) 

with Ji(n, n') and J2(n,n') given in Appendix [BJ 

Since the work done by the force F^^\t) (on the mov- 
ing wall) from a Fermi particle is supplied by the excess 
energy induced by its energy loss, we find 

The force can also be calculated by taking the expec- 
tation of the force operator F in Eq. ([22l) defined in the 
original space: 



/■^(*) ^ 1 ht 

Fit) = J ^*ix,t)F^ix,t)dx = —h + -^Im(/i) 



(71) 



Eqs. ([70)1 and (ITT]) exactly agree mutually, which guaran- 
tees the validity of the force operator defined in Section 

im 

To investigate the nature of the nonadiabatic force 
more carefully, however, we must estimate the integrals 
Iq and Ii . By expanding the exponential in Eq. ()65p as 

exp UHiQ)mv') =1 - f i(0)i(0)y^ 



-iLiQ)m)^V^ 



(72) 



we find (see Appendix [B|) 



7r2?i2(i(o)i(0))2 



In Ea. ((74|) the last factor can be taken as cos(- • • ) '-^ 1 
in the time range in Eg. (j49|) . We find that the i(0)- 
dependent terms are included not only in Ii but also in 
/q. Substituting Eqs. ^ and ^ into Eq. dTT]), we see: 
(i) the L(0)-independent term in Eq. ([7T|) gives rise to 
the adiabatic force iFad)] (ii) the remaining terms give 
the nonadiabatic force iFnon-ad)- 

Picking up the first term on r.h.s of Eq. (j73p . multiply- 
ing statistical weight /„ and summing up over all initial 
eigenstates, we find 



frij 



(75) 



which justifies Eq. 

Taking together Eq. ((7^ and the second term on r.h.s. 
of Eq. (j73p . multiplying statistical weight /„ and sum- 
ming up over all initial eigenstates, we find 



Fnon — ad 



c 



Ami 
m 



(76) 



with 



E ■ 



16^2 



(to2 



ivn? 



h^TT^ /^2 



2n27r4 



fn 



(77) 



Y,nAMnJ))^-i^J3ie,e) 

\ n 

(73) 



and 



Im(/i) =- ^i(0)L(0)^2nJi(n,^)J2(n,£) 



16 \ inlf 



7r2 ; (n2 - P-f 



X cos \ —il- n^)r(i) 



(74) 



which is again negative due to the dominant term propor- 
tional to V? . Irrespective of the direction of the moving 
wall, the non-adiabatic force always acts inwards and is 
proportional to the square of the wall velocity, which is 
in marked contrast with the general belief that the non- 
adiabatic force should be linear in the wall velocity and 
mimic the irreversible or frictional force. There is a mi- 
nor discrepancy between the absolute values of C and C", 
which is due to the difference in the way of solving the 
problem 23 1. 

\ In closing this Section, we should note the following 
iwo remarks: 

i) Firstly there is a mystery in obtaining the nonadi- 
abatic force in the exact analysis above. In this Section 
we had recourse to the transitionless states as basis func- 
tions. We can see neither nonadiabatic transition nor 
nonadiabatic force so long as tracking individual transi- 
tionless states. In fact, if we shall evaluate the expecta- 
tion value F = using only a single transitonless state 
in Eq. ()61|) . we will obtain formally the same result as 
in Eq. ([711) but with /q = 7r2£2 and Im(/i) = and can 
see no nonadiabatic force. Throughout this paper, we are 
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considering the case that the wall is fixed up to the initial 
time t = and suddenly moves at t > 0. Therefore the 
eigenstate under the fixed boundary generates at i = a 
mixture of the transitionless states that are eigenstates 
of the moving boundary, giving rise to nonvanishing co- 
efficients {c„}. Exploring Eqs. ([67|) and ([68]) . we can un- 
derstand that the correlation among non-zero coefficients 
{c„} resulted in the nonvanishing nonadiabatic force in 
Eq. ([76|) . By contrast, it is quite easy to see the mecha- 
nism for nonadiabatic force in the perturbative theory of 
Section IIIIl where the energy diffusion among standard 
adiabatic states can explain the nonadiabatic force. 

ii) Secondly the exact analysis can also reveal the na- 
ture of the non-adiabatic force in the case when the ex- 
panding rate of the cavity is not constant, so long as 
L{t) obeys Eg. ((5^ . namely when L{t) = V at^ + bt + c. 
Under the initial condition with L(0) = Lq and L{0) = 
2^ 7^ 0, F{t) can be calculated, leading to the identical 



result Fn 



C 



-LloT 



with a negative constant C. 



The details are given in Appendix |3 Thus our assertion 
that the non-adiabatic force is quadratic in the wall ve- 
locity and thereby does not break the time-reversal sym- 
metry does hold also for the hard-wall cavity with the 
time-dependent wall velocity. 



V. CASE OF SOFT- WALLED CONFINEMENT 

To see the universality of our argument so far, we pro- 
ceed to investigate the case of the soft- wall confinement, 
and consider the force acting on the soft wall. Here the 
Fermi gas is assumed to be confined in a harmonic trap 
with the confining length L changing linearly in time (see 
Fig©. 

To evaluate the expectation of energy and force, we 
shall solve the Schrodinger equation for a particle un- 
der the harmonic trap with the time-dependent trapping 
frequency, 



ih 



dt 



Hip{x,t) 



2 dx^ 



+ -muj {t)x ip{x,t), 



-OO < X < +00, 



(78) 



where, due to the prescription h^/m ~ 1, we see m = 
h^. In Eq. ([75| uj{t) is expressed in terms of the time- 
dependent confining length: huj{t) — -jT(t) where L{t) = 
Lq + Lt with constant L. 

After scale and gauge transformations like Eqs. ([55)) 



■L(t)- 



\ 



X 



FIG. 2: Soft wall confining the quantal gas. L{t) is the time- 
dependent confining length. 



and dST]), Eq. ^ is reduced to 



ihL(t) 



dt 



2 dy^ 



+ -h'y'^{y,t), (79) 



where ip{x,t) — ,^ ex.p{^hL{t)L{t)y'^)(l){y,t) and x 
L{t)y. The solution of Eq. ([79]) can be written as 



(j)„{y,t) = exp -i 



ftL2(4 



ds Y,i{y), 



(80) 



where 



r„(,) = J^exp(-^)i74V)^y) (81) 



with 



i7„(z) = (-l)"e^ — e- 



(82) 



An assembly of solutions in Eq. (|80|) . which are transi- 
tionless states, constitute the complete ortho-normal set. 

As in the case of the hard- wall cavity, let's require that 
the initial state is an eigenstate under the fixed harmonic 
trap: 



1 t X 

'ip{x,t)\t=o =ipoix) = . Yf 



Vi(0) 



(83) 



which is equivalent to 



0(y, 0) = Yi{y) exp i^--hL{0)L{0)y^ j . (84) 

The time-dependent solution satisfying this initial condi- 
tion is expressed as 



-t-oo 



Hy^t) = '^Cn<l>n{y,t), 



(85) 
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Expansion coefficients c„ are determined by 

Xj 

^{y,Q)Y:{y)dy 



— OO 
+ 00 



r,(y)r„(y)exp(--ni(0)i(0)y2)dy. (86) 



The average energy {E{t)) is calculated by 

/+00 
-00 

= ^K, + ^lMK,) + ^-K,, (87) 

where K, ^ /+^ (I'^.P + ^1 = /+r#*</),dy 

andi^2^/+^y2|^|2^y. 

The average force is obtained by taking the derivative 
of {E{t)) with respect to L as 



On the other hand, one should evaluate the expecta- 
tion value of the force operator to reproduce the above 
result. Although the range of x is not limited in the case 
of soft-wall confinement, one can define the force opera- 
tor using the time-dependent canonical transformations 
related to scaling of both coordinates and wave functions 
and its inverse transformations in Section |lll The force 
operators in the original space is given by 



F = ~- 



1 



hL 



d 



1 



(89) 



which includes the second term missing in Ea.(|22p. 
Therefore its expectation is 



l^Ko + ^lia{Ki), (90) 



F={F) 



which accords with Eq. and guarantees the validity 
the definition of the force operator in Eq. (|89|) . 

With use of asymptotic expressions for A'o and Ki, 
we again find Fnon-ad = C^o/t ^l{o) with a coefficient 
Csoft < , namely, the nonadiabatic force is proportional 
to the square of velocity of the soft wall, and never breaks 
the time-reversal symmetry. 



VI. CONCLUSIONS 

We investigated the nature of the non-adiabatic force 
acting on the cavity wall, which is generated from the 
non-interacting quantal gas confined in the expanding 



cavity. Firstly, with use of the time-dependent canon- 
ical transformations by which we can move to the non- 
expanding cavity, the force operator is defined. Secondly, 
we analyzed the expectation of the force operator per- 
turbatively with use of von Neumann equation for the 
density operator, which works when the cavity wall sud- 
denly begins to move at time origin. We found that the 
non-adiabatic force is quadratic in the wall velocity and 
thereby does not break the time-reversal symmetry, in 
marked contrast with the existing conjecture. Finally, 
using an assembly of the transitionless quantum states, 
we obtain the nonadiabatic force exactly. The exact re- 
sult justifies the validity of both the definition of force 
operator and the issue of the perturbative theory, and 
guarantees the present findings in the general case when 
the expansion rate of the cavity is time-dependent. The 
mysterious mechanism of nonadiabatic transition with 
use of transitionless quantum states is also explained. 
The study is done on both cases of the hard-wall and 
soft-wall confinement with the time-dependent confining 
length. Quantum fiuctuation theorem, deviation from 
the standard Fermi-Dirac distribution and equation of 
states, etc. in the expanding cavity where Hilbert space 
is time-dependent also constitute interesting subjects, 
which will be investigated in due course. 
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Appendix A: The case L{t) = \/ at^ + bt + c 

We consider the expanding cavity of the size L{t) gov- 
erned by Eq. ([59|) . which has the general solution L{t) = 
\J at'^ +bt + c with = — 4ac and the initial velocity 
L(0) = ^ ^ 0. In this case the reduced Schrodinger 



(Al) 



equation in Eq. ([58|) takes the following form: 

..j-2d(t) ld^(j) 1 2p2 2, 

inL —— — -—r h B y (h, 

dt 2dy^ 8 " ^' 

where y = x/L{t) and 

0(y, t) = VIexp (^-^L V) Hvm.t). (A2) 
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Changing the time variable from t to r defined by 
Eq.dnH), Eg. ([XT]) can be solved as 



= ¥'„(2/)exp ( -KnT{t) 



(A3) 



with 



(Pn{y) = An exp ( — I y Re 



, ^ 1 3 Kn 3 ihB , 
' 4 fiS' 2' 2 ^ 
(A4) 



where M{a,b,y) is Kummer function (i.e., the confluent 
hypergeometric function) and if„ are the solutions of the 
equation 



^.^ /3 3 z/iB 

iVi hi — , — . — 

A hB' 2 2 



0. 



(A5) 



In the semiclassical region where hB ~ 0, one can see 

2 2 

Now we shall solve the time-dependent problem under 
the initial state 



</)(y,0) = V2sin(^£y), 



(A6) 



2 sin(f^) with the 



which corresponds to ipix, 0) — y -^py olln^-^py 

eigen-energy EAO) = ^^ttt- Expanding (jjiy, t) in terms of 
the complete set of functions (/?„(y), we have 



where 



<l>{y,'t) = ^Cn{t)Lpn{y), 



Cn{t) = C„(0)exp ( --KnT{t) 



with c„(0) given by 
1 

c„(0) = x/2 j ^;(y) sinfTry 



=Jo(n,€)-^L(0)L(0)J2(n,^) 

-^(L(0)i(0))2j3(n,£) + ... . 

Jo, J2 and J3 are defined respectively as 
1 

Jo(n,^) = V2 y <(y) sm{n£y)dy, 




(A7) 



(A8) 



exp 1 -^L(0)L(0)2/^ ) dy 



(A9) 



(AlO) 



1 

J2inJ) = V2 j yV:(2/)sin(^£y)dy, (All) 



and 



J3(n, £)^V2j y^ifiUy) sin(7r£y)dy. (A12) 


Pair products of c„ s are 

c„-(0)<(0)- Jo(n',£)Jo(n,£) 
in • 

-yL(0)L(0) [Jo(n',£)J2(n,f) - Mn,£)J2{n',i)] 

-^(L(0)L(0))2 [Jo(n',^)J3(n,^) 
+Jo(n,^)J3(n',^) + 2J2(n',^)J2(n,^)] (A13) 

and 



Cn'it)cl{t) = c„'(0)c*(0)exp ~i 



-ds 



(A14) 



The expectation for the force operator in Eg . ([221) is 
obtained as 



L(t) 



(A15) 



(F^^^(t)) = / il;*{x,t)F^{x,t)dx 
Jo 

= Jl^o + -^Im(/i), 

with Jo and Ji being expressed as 
1 

Io = J |0„pdy = 2^|c„(O)|2i^„ 


1 

+ -h^B^ I] I] Re [c:(t)c„.(i)] / Vn^wy^dy 

(A16) 



and 



1 

/ 



h = I y(p*(j)ydy 



= EE<(*K'W/ y^^dy, (A17) 

respectively. Using the asymptotic forms for the pair 
products in Eg. (|A13p together with Eg. (|A14p and taking 
the short-time region employed in the main text, we reach 
the result: 



Fno7i—-ad — C 



m ' 



(AI8) 



which is again proportional to the square of the wall ve- 
locity. 



12 



Appendix B: Calculation of / and J 



Ji , J2 and J3 used above are given respectively by 



Coefficients c'n (0) are defined by 



(0) 



=2 / sin (/Try) sin {mry) exp 
Jo 

=<5„£- yi.(0)L(0)J2(n,^) 
_^(i(0)i(0))2j3(n,£) + 
Their products are 

ci'?*(0)cW(0)^<5„-,<5„, 



-^L(0)L(0)y2)dy 



(Bl) 



Ji{n,£) — 2 / y sm{£Try) cos{mry)dy 
Jo 



(-1) 



TT^, for n = 



otherwise, 



J2(n,£) = 2 / y sin(£7ry) sin(n7ry)c?y 
"'0 

( 11 



for n ~ £ 
(-l)"+ \„._%.^. otherwise 



(B4) 



(B5) 



- yi(0)L(0) [(5„,,J2(n,€) -5„,J2(n'/)]- 
- :^(i(0)i(0))2 [Sn,,Mn,£)+6^iMn',£) 
+2J2{n',£)J2{n,£)]. (B2) 

Therefore /i is expressed as 



J2 nJiin, £)J2(.n, £)-£j2 n)J2{n', £) 



(B3) 



and 



J3{n,£) = 2 / y sin(£7ry) sin(n7ry)c?y 
Jo 



5 2n*7r-* ' 

for n = £ 



otherwise. 



(B6) 
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